Relativistic descriptions of final-state interactions in charged-current quasielastic 
neutrino-nucleus scattering at MiniBooNE kinematics. 
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The results of two relativistic models with different descriptions of the final-state interactions are 
compared with the MiniBooNE data of charged-current quasielastic cross sections. The relativistic 
mean field model uses the same potential for the bound and ejected nucleon wave functions. In the 
relativistic Green's function (RGF) model the final-state interactions are described in the inclusive 
scattering consistently with the exclusive scattering using the same complex optical potential. The 
RGF results describe the experimental data for total cross-sections without the need to modify the 
nucleon axial mass. 
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The double differential cross sections for muon neu- 
trino charged-current quasielastic (CCQE) scattering, re- 
cently measured by the MiniBooNE collaboration [l[, 
have raised debate over the role of the various theoret- 
ical ingredients entering the description of the reaction. 
High-quality descriptions of the CCQE differential cross 
sections in the few-GeV region are required to support 
neutrino oscillation measurements [l|, y| . The energy re- 
gion explored requires a relativistic description of the pro- 
cess, where not only relativistic kinematics is considered, 
but also nuclear dynamics and current operators should 
be described within a relativistic framework. 

The simplest relativistic model to describe CCQE neu- 
trino scattering is the relativistic Fermi Gas (REG). 
When a dipole shape is assumed for the axial form factor, 
the nucleon axial mass Ma has been used as a free pa- 
rameter within the RFG model. Indeed, the MiniBooNE 
cross section [1| is underestimated by the REG unless 
Ma is significantly enlarged (1.35 GeV/c^) with respect 
to the accepted world average value (1.03 GeV/c^ [3|]). 
As it turns out from comparison with electron scattering 
data, the REG is too crude to correctly account for the 
nuclear dynamics. Thus, a larger axial mass within the 
REG could be a way to effectively incorporate nuclear 
effects. 

More sophisticated models have been applied to 
neutrino-nucleus scattering. At the level of the impulse 
approximation (lA), models based on a realistic spectral 
function [^ or on the relativistic lA (RIA) which, con- 
trarily to RFG, are in good agreement with electron scat- 
tering data, also underestimate the experimental CCQE 
cross sections p-uf unless Ma is significantly enlarged, 
as indicated by the REG prediction. 

It has been pointed out that in some kinematic re- 
gions where the neutrino flux for the experiment has 



significant strength, the reaction may have sizable con- 
tributions from effects beyond the I A. These include 
two-particle-two-hole (2p-2h) excitations, which may be 
reached via two-body meson-exchange currents (MEG). 
The contribution of the vector MEC in the 2p-2h sec- 
tor, evaluated in the model of Ref. [8[, has been incor- 
porated in a phenomenological approach, indicated as 
SuSA, based on the SuperScaling behavior of electron 
scattering data [y]. The strict SuSA predictions show 
a systematic discrepancy in comparison with the Mini- 
BooNE cross sections [g|- The inclusion of 2p-2h MEC 
contributions yields somewhat better agreement with the 
data [71] , although theory still lies below the data at larger 
angles. Other theoretical results that incorporate multi- 
ple knockout excitations are in accordance with the ex- 
perimental cross sections without the need to increase 
the value of Ma id, llll ■ We note that the results in [10| 
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are obtained in a relativistic approach and in 
non-relativistic one. 



Within the QE kinematic domain, the treatment of the 
final-state interactions (ESI) between the ejected nucleon 
and the residual nucleus has been proved to be essential 
to compare to data. The relevance of ESI has been clearly 
stated for exclusive {e,e'N) processes, where the use of 
complex optical potentials in the distorted-wave impulse 
approximation (DWIA) is required |12h14|. In the anal- 
ysis of inclusive reactions, FSI remains a crucial ingredi- 
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ent for a proper description of data 
and inelastic channels contribute to the inclusive process. 
Thus, the complex potential, with imaginary terms de- 
signed to reproduce just the elastic channel, should be 
dismissed. Different approaches have been used to ac- 
count for FSI under inclusive conditions. For instance, 
in the approaches based on the relativistic DWIA (RD- 
WIA), they have been accounted for by using purely real 



potentials. The final nucleon state has been evaluated 
with the real part of the relativistic energy-dependent 
optical potential (rROP), or with the same relativistic 
mean field potential considered in describing the initial 
nucleon state (RMF) [11, [l^. However, the rROP is un- 
satisfactory from a theoretical point of view, since it is an 
energy-dependent potential, reflecting the different con- 
tribution of open inelastic channels for each energy, and 
under such conditions dispersion relations dictate that 
the potential should have a nonzero imaginary term J22l | . 
On the other hand, the RMF model is based on the use 
of the same strong energy-independent real potential for 
both bound and scattering states. It fulfills the disper- 
sion relations 2^ and also the continuity equation. The 
RMF model applied to inclusive QE (e, e') processes de- 
scribes scaling behavior and gives rise to a superscaling 
function with a sig nificant asymmetry, in good agreement 
with data [H, [ij 

A different description of FSI makes use of relativis- 
tic Green's function techniques 
suitable approximations [if 
cally related to the lA, the components of the nuclear 
response are written in terms of the single-particle (s.p.) 
optical model Green's function. This formalism allows 
to recover the contribution of non-elastic channels in the 
case of inclusive scattering, starting from the complex rel- 
ativistic optical potential (ROP) which describes elastic 
nucleon-nuclcus scattering data. The relativistic Green's 
function (RGF) model allows for a consistent treatment 
of FSI in the exclusive and in the inclusive scattering and 
gives also a good description of (e, e') data 16l.ll7|. 

The results of the RMF and RGF models have been 
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compared for the inclusive QE electron scattering [17 1 
and for the CCQE neutrino scattering [27[. As men- 
tioned, both models describe successfully the behavior 
of electron scattering data and their scaling and super- 
scaling functions and both produce a significant asymme- 
try in the scaling function that is strongly supported by 
data. There are, however, some differences between the 
RMF and RGF results depending on kinematics, which 
increase with the momentum transfer. Whereas the RMF 
may be considered as a faithful representation of the pure 
"nucleonic" contribution to the inclusive response, the 
RGF, on the contrary, may to some extent translate loss 
of elastic strength to non-nucleonic degrees of freedom, 
contributing to the imaginary optical potential, into in- 
clusive strength predicted by the RGF. 

In this letter the predictions of the RMF and RGF 
models are compared with the recent CCQE MiniBooNE 
data. The comparison between the results of the two 
models [l7|, [27 [ , which make use of very different ingredi- 
ents, can be helpful for a deeper understanding of nuclear 
effects, more specifically FSI, which may play a crucial 
role in the analysis of CGQE data and its infiuence in 
studies of neutrino oscillations at intermediate to high 
energies. This is of particular interest for the case of the 



MiniBoone CCQE data which, given the nature of the 
experiment, rnay receive more than pure nucleonic con- 
tributions [28[. Thus the RMF would represent a lower 
bound to MiniBoone CCQE, while the RGF should yield 
larger predictions. 

Details of the two models can be found in ll{ 
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|29|, 130| for the RMF and in [li 
the RGF. In the RMF case, the components of the nu- 
clear response are obtained from the sum over all the 
s.p. shell-model states of the squared absolute value of 
the transition matrix elements of the single- nucleon cur- 
rent. In the RGF case, the calculations require matrix 
elements involving the eigenfunctions of a complex opti- 
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cal potential and of its Hermitian conjugate [If 

In both calculations the bound nucleon states are self- 
consistent Dirac-Hartree solutions derived within a RMF 
approach using a Lagrangian containing a, w, and p 
mesons [33|. The same real potential gives the scat- 
tering states in the RMF, whereas in the RGF calcu- 
lations two parameterizations for the ROP have been 
used: the energy-dependent and A-dependent EDADl 
and the energy-dependent but A-independent EDAI-12C 
complex phenomenological potentials of |34{ . which are 
fitted to proton elastic scattering data on several nuclei 
in an energy range up to 1040 MeV. The comparison be- 
tween the results obtained with two different phenomeno- 
logical optical potentials may indicate how the incom- 
plete determination of this important ingredient can in- 
fiuence the predictions of the model. In all the calcu- 
lations we have used the standard value of the nucleon 
axial mass, i.e., Ma = 1.03 GeV/c^. 

In Fig. [U we show the CCQE double-differential ^^C 
(i'^, /i~) cross section averaged over the neutrino flux as 
a function of the muon kinetic energy T^ . In each panel 
the results have been averaged over the corresponding 
angular bin of cos 6, where 9 is the scattering angle of the 
muon. The results evaluated with RMF (blue line) and 
RGF with EDADl (red) and EDAI (green) potentials are 
compared with the MiniBooNE CCQE data [l[. 

The RMF results [7| yield reasonable agreement with 
data for small angles and low muon energies, the discrep- 
ancy becoming larger as 9 and T^ increase. The shape 
followed by the RMF cross sections fits well the slope 
shown by the data. A good agreement with the exper- 
imental shape is shown also by the RGF cross sections. 
The RMF and RGF models yield close predictions at 
larger values of T^ for all the bins of cos 9 shown in the 
figure. Notice, however, that the RGF cross sections are 
generally larger than the RMF ones, particularly around 
the peak region, where the RGF produces cross sections 
in reasonable agreement with data. 

It is worth noticing that the differences between the 
RGF results obtained with the two optical potentials are 
enhanced in the peak region and are in general of the 
order of the experimental errors. The EDADl and EDAI 
potentials yield close predictions for the bin 0.4 < cos 9 < 
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FIG. 1. (color online) Flux-averaged double differential cross 
section per target nucleon for the CCQE ^^C(!^,j,/i^) reac- 
tion calculated in the RMF (blue line) and in the RGF with 
EDADl (red) and EDAI (green) potentials and displayed ver- 
sus T^ for various bins of cos 9. The data are from Mini- 
BooNE [l|]. The uncertainties do not include the overall nor- 
malization error 5A'^=10.7%. 



0.5, a small differences is seen in the bin 0.8 < cos 9 < 0.9, 
being the RGF-EDAI cross section larger than the RGF- 
EDADl one, while the difference is sizeable for the bin 
0.7 < cos6' < 0.8, with the RGF-EDADl resuhs closer 
to the RMF than to the RGF-EDAI ones. 

The RMF model uses the effective mean field, that re- 
produces the saturation behavior of nuclear matter and 
the properties of the ground state of nuclei. It includes 
only nucleonic contributions to the inclusive process. The 
RGF uses phenomenological optical potentials, fitted to 
elastic proton-nucleus scattering. The loss of elastic flux 
into inelastic channels (either multi-nucleon knockout as 
well as non nucleonic excitations) caused by the imagi- 
nary term of these potentials is recovered for the inclu- 
sive scattering making use of dispersion relations. The 
larger cross section shown by the RGF can be attributed 
to non purely nucleonic inelasticities represented in the 
phenomenological ROP [H, ^ . 

In Fig.[2]the flux-averaged double differential cross sec- 
tions are plotted versus cos 9 for two bins of T^, i.e., 
0.2 < Tf, < 0.3 GeV and 0.6 < T^ < 0.7 GcV. The 
approximate shape of the experimental cross section is 
well described by the models. The RMF results gen- 
erally underestimate the data, especially for the lower 
muon energy values, the agreement improves as T^ in- 
creases. The RGF provides a better accordance with the 
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FIG. 2. (color online) Flux-averaged double differential cross 
section per target nucleon for the CCQE ^^C{u^, fj,~) reaction 
displayed versus cos 6 for two bins of T^. The results obtained 
with RMF (blue line), RCF EDADl (red), and RGF EDAI 
(green) potentials are compared with the MiniBooNE data 
of[l. 



size of the experimental cross section. The agreement is 
better for smaller angles while the data are slightly un- 
derpredicted as 9 increases. The RGF-EDADl yields in 
general a lower cross section than the RGF-EDAI, yet 
higher than the RMF one. 
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FIG. 3. (color online) Total CCQE cross section per neutron 
versus the neutrino energy. The cross sections calculated in 
the RMF (blue line), RGF EDADl (red), and RGF EDAI 
(green) potentials are compared with the flux unfolded Mini- 
BooNE data of 0. . 

Finally, in Fig. [3] the total QE cross section per neu- 
tron obtained in the RMF and RGF models are displayed 
as a function of the neutrino energy Ei, and compared 
with the "unfolded" experimental data [1| . It was shown 
in [7[ that the differences between the results of the RMF, 
SuSA, and rROP models tend to be washed out in the 
integration and that all these models, representing es- 



scntially the same nucleonic contribution to the inclu- 
sive cross sections, yield very similar results, all of them 
underpredicting the total MiniBoone CCQE experimen- 
tal cross section. Larger cross sections, in particular for 
larger values of E,^, are obtained in the RGF with both 
optical potentials. The differences between RGF-EDAI 
and RGF-EDADl are here clearly visible, being RGF- 
EDAI in good agreement with the shape and magnitude 
of the experimental cross section and RGF-EDADl above 
RMF but clearly below the data. The differences between 
EDAI and EDADl are due to the different values of the 
imaginary parts of both potentials, particularly for the 
energies considered in kinematics with the lowest 6 and 
the largest T^. These kinematics, which were not con- 
sidered in previous RGF calculations, give large contri- 
butions to the total cross section and emphasize the dif- 
ferences between the RGF predictions with both optical 
potentials. Notice that EDAI is a single- nucleus parame- 
terization, which docs have an edge in terms of better re- 
production of the elastic proton- ^^C phenomenology |34| 
compared to EDADl. and also leads to CCQE results in 
better agreement with data. 

Summarizing, in this letter the results of the RMF and 
RGF models have been compared with the recent CCQE 
MiniBooNE data. Both models give a good description 
of the shape of the experimental cross sections. The 
RMF generally underpredicts the data, particularly for 
lower values of 9 and T^. In contrast, the RGF can give 
cross sections of the same magnitude as the experimen- 
tal ones without the need to increase the standard value 
of the axial mass. The larger cross sections in the RGF 
model arise from the translation to the inclusive strength 
of the overall effect of inelastic channels (nucleonic and 
non- nucleonic) . At present, lacking a phcnomenological 
optical potential which exactly fullfills the dispersion re- 
lations in the whole energy region of interest, the RGF 
prediction is not univocally determined from the elas- 
tic phenomenology, though some preference to the EDAI 
predictions should be given. 

Our results give a further and clear indication that 
before drawing conclusions about the comparison with 
data, a careful evaluation of all nuclear effects and the 
possible effect of some non-nuclenic contributions to 
CCQE MiniBoone data is required [28[. This is impor- 
tant also to reconcile former results for which RMF was in 
good agreement with previous CCQE data [35[ A better 
determination of a phcnomenological relativistic optical 
potential which closely fullfills the dispersion relations 
deserves further investigation. 
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